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We introduce the map representation of a time-delayed system in the presence of delay time 
modulation. Based on this representation, we find the method by which to analyze the stability of 
that kind of a system. We apply this method to a coupled chaotic system and discuss the results in 
comparison to the system with a fixed delay time. 
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In a real situation, time delay is inevitable, since the 
propagation speed of an information signal is finite. The 
effects of time delay have been investigated in various 
fields in a form of delay differential equation 0, 0, EJ Dl 
5, 6]: 



x = F(x(t),x(t - ?)), 



(1) 



where f is the delay time. Since the Volterra's predator- 
prey model [ij, delay time has been considered such 
various forms as fixed HHHH, distributed^], state- 
dependent [j| , and time-dependent 0, 0, ones. Up to 
now, the delay effects of these forms in dynamical systems 
have been extensively studied not only in various helds 
of physics 0, biology 0, and economy p| but also for 
the cases of application^, El 0] ■ The recent discovery 
@ that the neural synchrony is enhanced by the discrete 
time delays shows the delay feedbacks play a significant 
role in synchronization phenomenon. This discovery calls 
new scientific attention on the time-delay feedbacks and 
synchronization phenomenon. 

Meanwhile, in regard to the application of time- 
delayed systems to communication, delay time modu- 
lation (DTM) H IH |H El driven by the signal of a 
independent system was introduced to erase the imprint 
of the fixed delay time. It was reported that the DTM 
increases the complexity of the time-delayed system as 
well as prevents the system from being collapsed into the 
simple manifold in reconstructed phase space, which was 
a serious drawback of the time-delayed system to the ap- 
plication to communication. It was also found that the 
synchronization in the presence of DTM can be estab- 
lished and the death phenomenon is enhanced, due to 

dtm[S Sl- 
it was further shown that the time-delayed system can 
be represented by high dimensional coupled maps 
and the methods by which to understand the charac- 
teristics of the time delayed system with discrete time 
delays have been developed ■ Accordingly, the natural 
questions to be tackled have become " how can we find the 
coupled-map representation of the time-delayed system in 
the presence of DTM ? " and " how can we analyze the 
stability of that kind of system T . The answers on these 



two questions are quite essential to understand the syn- 
chronization and the enhancement of death phenomenon 
mentioned above with regard to DTM |HJ,[l3|]. In this pa- 
per, we newly introduce the coupled-map representation 
of a time-delayed system in the presence of DTM, and 
develop a method which enables us to get the Lyapunov 
exponents of the system based on the representation. We 
shall also show the DTM accelerates the transition to hy- 
perchaos in a coupled time-delayed system by presenting 
the Lyapunov spectrum obtained by the method. 

The delay differential equation of Eq. (1) can be 
rewritten by the discrete form: 



(2) 



where x n is the value of x at n-th time step such 
that x n = x(nAt) and x n - T — x(nAt — tA£). Here 
f(x n ,x n - T ) = x(nAt) + F(x(nAt),x(nAt - rAf))At. It 
is known that the above system is equivalent to the (t+1) 
dimensional coupled maps [T^ : 



X n+1 




X 7l+1 


- x T - 1 




- x T ~ 2 



(3) 



"n+l 



where the variables a; 1 , . . . , x T n construct the "echo type" 
feedback loop. The above representation plays an impor- 
tant role in analyzing the time-delayed system with fixed 
delay times. DTM jj consists of a modulation system, 
y = G(y) which generates the driving signal y„, and a 
scaling system, f = H(y) which adjusts the amplitude 
of modulation. Accordingly two additional equations are 
introduced: 



Vn+l 



g{y n ), 



(4) 



where g(y n ) = y(nAt) + G(y(nAt)) and h(y n ) = 
H(y(nAt)). The modulated delay time r in Eq. (4) 



2 



selects the length of delay time for a feedback variable 
x„_ r . As shown in Eq. (3) the feedback variable (say 
x^) is not changed in the fixed delayed system. However, 
in the time-dependent delayed system, we can interpret 
that the feeding variable is selected according to the delay 
time (say x\ with k — 0, 1,2, ...,r m ) in Eq. (3). There- 
fore we introduce the map representation as follows: 



v n+l 



"n+1 



f(x n7 X Tl 
™T™-1 



(5) 



k 

Vn+l 



HVn), 

g{yn), 



(6) 



where r m is the maximum of the amplitude of DTM and 
k is an integer number which determines the variable to 
be fed into the systems. We have R Tm+1 x I x R 1 phase 
space. From Eq. (5) and (6), we get the Jacobian matrix 
as follows: 



df(x n x n ) df(x° n .x k n ) cr m df(x° n ,x n ) cr m -l df(xl,x*) <-T„ 
dxl dxl °k dx>* °k dx>* °k 



°k 



dk 











where Si = 1 only if k = I else Si = 0. We introduce 
the (r m + 3) dimensional orthonormalizcd initial vectors 
such that: yi = (1, 0, 0, . . . ), y 2 = (0, 1,0,...), and . . . 
which are evolved by y^ = Jy^ . By following the standard 
procedure |l6j , we can find the Lyapunov spectrum of the 
time-delayed system in the presence of DTM as follows: 

A <=^i>^> (8) 

fc=l 

where r is the number of orthonormalization within the 
chosen finite time interval T and is the norm of the 
i-th vector at fc-th orthonormalization process. 

As an example of DTM, we consider the two coupled 
logistic maps as follows: 

X n +1 — lX n (\ — X n ), 

t = [Ay n ], 

Vn+l = 70?/n(l - Vn), (9) 

Here x r n — (1 — a)x n + a{xn-r) and a is a coupling 



strength. We took 70 = 4.0 and A = 4.0 and [x] denotes 
the maximum of the integer which is less than x and r m = 
max([Aj/ n ]) = 3. According to the general representation 
of Eqs. (5) and (6), the system described in Eq. (9) 
can be rewritten by the 6-dimensional coupled maps as 
follows: 



4+i = lx k n {l-x k n )=F{x k n ), 

X n+1 ~~ x n> 

x n+l = x ni (10) 

k = [Ay„], 

Vn+l = 702/n(l - Vn), (H) 

where x\ — (1 — a)x^ +ax k l . The corresponding Jacobian 
matrix is given by: 



3 



7 ((l-a) + a/3<$)(l 



1 






2x k n ) yad k (l-2x k ) jaS k (l-2x k ) 1 a5 k {l-2x k n ) F(x k+1 ) - F(x k ) 









7o(l - 2y n ) 



(12) 



where we have used the relation — gfc' = 
F(a:°,a;* +1 ) - F(af°,a;*) noting that k is an integer [l7| 
and Here d ^"^ is not zero only at 0.25, 0.50, and 0.75 
which arc unstable periodic orbits of the driving logistic 
map (the second equation of Eq. (11)). Therefore we put 
a ^y"^ = to the above equation. 

To verify our method, we calculate the Lyapunov spec- 
trum in the case of DTM as well as in the fixed delay 
case. We can easily restore Eq. (10) to the fixed delay 
case by fixing k = 3 without Eq. (11). One can see 
the bifurcation diagram presented in Fig. 1 (a) and the 
corresponding Lyapunov spectrum in Fig. 1 (c). The 
quasi-periodic regime in 7 6 [3.1,3.7] and the periodic 
window in 7 G [3.7, 3.8] are exactly matched with the 
corresponding Lyapunov spectra. As one can see, the bi- 
furcation diagram is different from that of the single lo- 
gistic map, which is due to the delay feedback |l5j. In the 
case of fixed delay (Fig. 1 (a) and (c)), two largest Lya- 
punov exponents become positive above 7 = 3.8 through 
the quasi-periodic regime in 7 6 [3.1, 3.7]. That is to say, 
the system develops to hyperchaos above that point . 

The Lyapunov exponents in the case of DTM are pre- 
sented in Fig. 1 (d). Actually we have two more Lya- 
punov exponents which are not plotted in the figure. One 
has a constant value, A5 = 0.691 . . . , which describes the 
dynamics of isolated master system and the other has a 
relatively large negative value, \§ € [—0.5,-0.6]. While 
in fixed delay time, x r n is directly fed into the delayed 
system, in DTM all delay variables 2° , x x nl . . . x T n can con- 
tribute to the system depending on the state of the mod- 
ulating signal y n . For this reason, the largest Lypaunov 
exponents are very closer with each other (Fig. 1 (c)) 
than those of the fixed delay time (Fig. 1 (d)). The sum 
of positive Lyapunov exponents X)a >o ^» * s ^ ne en t r0 Py 
of th e sy stem which quantifies the complexity of the sys- 
tem |l8j . For example, at 7 = 3.95 the sum of positive 



Lyapunov exponents is 0.262 in the case of DTM, while 
it is 0.128 in the case of fixed delay time. This type of 
behavior caused by DTM is fulfilled with the result of the 
recent report in which the entropy is increased depending 
on the property of the driving signal for DTM [jjj- 

Figure 2 (a) shows the time series at the intermittent 
regime (7 = 3.5) and (b) shows the time series at the 
hyperchaotic regime (7 = 3.95) when the delay is modu- 
lated. The return maps are presented in Fig. 2 (c) and 
(d), which show a complex structure. Comparing the cor- 
relation functions in Fig. 2 (e) (with the fixed delay time) 
and (f) (with DTM), one is confirmed that the DTM re- 
duces the length of correlation |10|. We emphasize that 
all the procedure can be applied to another system re- 
gardless of coupled maps or flows in principle because in 
case of flow, the system also can be represented by the 
same form of coupled maps of Eqs. (5) and (6) (one just 
needs to note that the dimension of the phase space is 
determined by the value of delay time r and the chosen 
time step At such that dim ~ t/ At). 

In conclusion, we have introduced the coupled-map 
representation of the time-delayed system in the presence 
of delay time modulation and found the method by which 
to analyze the stability of the system. We have applied 
this method to a coupled chaotic system and confirmed 
that the Lyapunov spectrum of the time-delayed system 
in the presence of DTM can be exactly determined which 
explains the bifurcation behavior. Furthermore, we have 
observed that DTM erases the preference of the feedback 
variable, by scanning all delay variables during the evolu- 
tion, and that for this reason it leads Lyapunov exponents 
to be close to each other. 
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FIG. 1: The bifurcation diagrams and the Lyapunov spectra with a — 0.7. (a) and (c) show the bifurcation diagram and 
Lyapunov spectra at fixed delay time (only using Eq. (10) with k — r m ), (b) and (d) at DTM (using Eq. (10) and (11)). We 
have taken r = 5000 orthonormalization process with the period of T = 70 for each point (see Eq. (8)). 
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FIG. 2: Temporal behaviors and retrun maps of delayed lo- 
gistic map in the presence of DTM (Eq. (10) and (11)). (a) 
Time series of x„ at 7 = 3.5 (b) at 7 = 3.95; (c) return map 
x ( l versus x\ (d) x° n versus x\; autocorrelation functions of 
function of iteration at 7 = 3.95 in cases of (e) fixed 
delay and (f) DTM. 
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